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Global nitrogen cycling is mainly mediated by the activity of microorganisms. Nitrogen 
cycle processes are mediated by functional groups of microorganisms that are affected 
by constantly changing environmental conditions and substrate availability. In this study, 
we investigated the temporal and spatial patterns of nitrifier and denitrifier communities 
in an intertidal wetland. Soil samples were collected over four distinct seasons from three 
locations with different vegetative cover. Multiple environmental factors and process rates 
were measured and analyzed together with the community size and composition profiles. 
We observed that the community size and composition of the nitrifiers and denitrifiers are 
affected significantly by seasonal factors, while vegetative cover affected the community 
composition. The seasonal impacts on the community size of ammonia oxidizing archaea 
(AOA) are much higher than that of ammonia oxidizing bacteria (AOB). The seasonal 
change was a more important indicator for AOA community composition patterns, while 
vegetation was more important for the AOB community patterns. The microbial process 
rates were correlated with both the community size and composition. 
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INTRODUCTION 

The intertidal wetland ecosystem is a type of coastal wet- 
land that is characterized by high primary production rates 
and intense remineralization in its sediments (Carling, 1982). 
Intertidal wetlands and estuarine ecosystems are considered the 
most important land-sea interaction areas in the world. The estu- 
arine intertidal wetland is a natural barrier in purifying territorial 
pollution and attenuating the riverine load to the sea and plays a 
key role in controlling marine eutrophication and global nitrogen 
recycling. Microbial communities in these ecosystems are usually 
highly abundant and diverse (Wilms et al., 2006). 

The Dongtan wetland is an estuarine intertidal wetland near 
the sea that is located at the eastern end of Chongming Island in 
the Yangtze River estuary. Chongming Island is the third largest 
island in China and the largest alluvial island in the world (Gan 
et al., 2009). Because of its geographic location, organic mat- 
ter contained in this wetland is carried by both the Yangtze 
River and the open sea. The Yangtze River carries ~4.68 x 10 8 
tons of sediment to this area annually, and more than half of 
the sediment accumulates in the estuary region, which con- 
tains a large amount of nitrates, ammonium, and organic mat- 
ter (Chen, 1996) that are primary or secondary substrates for 
nitrogen transformation processes, including nitrification and 
denitrification. 

Nitrification and denitrification are two principal microbial 
processes that determine the nitrogen loss from ecosystems, both 
of which are mainly driven by the activities of microbial com- 
munities. Nitrification, including ammonium oxidization and 



nitrite oxidization, was long believed to be accomplished by a 
small specific group of bacteria, until the existence of an archaeal 
ammonium oxidizer was identified about a decade ago (Konneke 
et al, 2005). Conversely, denitrification has been detected in 
many different microbial taxa, most of which are phylogenetically 
distantly related (Zumft, 1992). 

In estuaries, the physiological and metabolic features of the 
aforementioned functional microbial communities in terms of 
activity, abundance, and composition are constantly changing 
due to the unstable environmental parameters, such as the C/N 
ratio (Li et al, 2012; Abell et al, 2013), pH (Cao et al, 2011), 
temperature (Sahan and Muyzer, 2008; Zheng et al., 2013), and 
salinity (Mosier and Francis, 2008; Santoro et al., 2008; Jin et al., 

201 1) . Unlike most studies in pelagic or terrestrial systems where 
the AOA always outnumber the AOB, there are some reports 
that have shown higher numbers for AOB than AOA in estuaries 
(Magalhaes et al, 2009; Wankel et al, 2011; Zheng et al, 2014); 
however, there are also some reports of higher AOA than AOB 
contents in different estuaries (Beman and Francis, 2006; Li et al., 

2012) . The actual contribution of AOB and AOA to the ammo- 
nium oxidization has also been the subject of intensive debate. 
There are some reports that showed a correlation of nitrification 
with AOB (Magalhaes et al, 2009; Wankel et al, 2011), while 
other reports suggest a high contribution of AOA to nitrification 
(Caffrey et al., 2007; Jin et al, 2011; Zheng et al., 2014). The rel- 
ative importance, community structure and abundance of AOB 
and AOA are complex issues that depend on actual environmental 
conditions (Berhard and Bollmann, 2010). 
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Thus, to better understand the in situ performance and 
function of nitrogen cycle microorganisms, including ammonia 
oxidizers and denitrifiers, dedicated studies of the community 
size and composition in relation to the environmental condi- 
tions and process rates are needed. To this end, we investigated 
the ammonia-oxidizing and denitrifying communities in the 
Dongtan wetland ecosystem and examined how the composi- 
tion and size of these communities correlate with environmental 
variables, such as season, vegetation types, and corresponding 
microbial process rates. 

MATERIALS AND METHODS 

SITE DESCRIPTION AND SAMPLE COLLECTION 

The Dongtan wetland lies in the typical subtropical monsoon 
region of China (Figure 1). The climate in this area is mild and 
characterized by four distinct seasons. The average annual tem- 
perature is 15.3°C, the average summer temperature is 26°C 
and the average winter temperature is 3°C (Xiao et al, 2009). 
Vegetation in the sampled site consisted of more than 90 types of 
vascular plants. Two types were dominant in both the distribution 
area and total biomass: exotic Spartina alterniflora and Phragmites 
australis. The low tidal zone of the Dongtan wetland is mudflats, 
also known as tidal flats, which have no vegetation cover (bare 
flat). The middle and high tidal zones are salt marshes. The mid- 
dle tidal zone is infested with Spartina alterniflora, and Phragmites 
australis thrives in the high tidal zone (Zheng et al., 2013). 

The samples were collected on 16 July 2008 (summer), 17 
October 2008 (autumn), 15 January 2009 (winter), and 17 March 
2009 (spring). During this period, the minimum environmental 
temperature was 6.0°C in January, and the maximum was 34.1°C 
in July (Table 1). At each sampling, three vegetative types of soil 
were sampled: soil in a bare flat (BF), soil covered by Spartina 



alterniflora (SA), and soil covered by Phragmites australis (PA) 
(Table 1). For each type of soil, four samples were collected with 
a cutting ring. Samples were then stored in a cooling box and 
immediately transported to the laboratory for further analysis. 
From each sample, 5-20 g of soil was used for measuring the 
soil's physico chemical properties, and 1-2 g was taken for DNA 
extraction. 

MEASUREMENT OF THE SOIL PROPERTIES AND MICROBIAL PROCESS 
RATES 

Soil pH and salinities were measured by the HANNA HI 9025 
pH/mV7°C waterproof meter (Hanna Instruments, Italy) using 
an extract of 1:10 dilution of soil with water. The concentrations 
of nitrate (NO^-N) and ammonium (NH^-N) were determined 
according to previously described methods (Solorzano, 1969; 
Norman et al., 1985). Total carbon (TC) and nitrogen (TN) 
were determined with a carbon, nitrogen, and sulfur analyzer 
(CNS-2000, USA) after combustion of the soil samples at 1250°C. 

Net nitrification and the potential denitrification rates were 
used as the microbial process rates in this study. The net nitrifica- 
tion rate was calculated based on the net production of NO^-N 
after incubating fresh soil samples (adjusted to 60% of the water- 
holding capacity) for 2 weeks aerobically in the dark at 25° C 
(Verchota et al, 2001; Rutigliano et al, 2009). 

The soil denitrification rate was determined as previous 
described (Yeomans et al, 1992). First, 20 g of fresh soil sam- 
ples were sealed in glass flasks with butyl rubber stoppers. Ten 
percent of the atmosphere contained in flasks was then replaced 
by C2H2to prevent reduction of N2O. After incubation at 25°C 
for 48 h, the potential denitrification rates were calculated by 
measurement of the N2O concentration using a gas chromato- 
graph (GC-14B, Shimadzu, Japan). 




FIGURE 1 | Map of the study area with sampling sites in the Dongtan wetland, eastern Chongming Island, China. BF, soil at bare flat (no vegetative 
cover); SA, soil covered by Spartina alterniflora; PA, soil covered by Phragmites australis. 



Frontiers in Microbiology | Terrestrial Microbiology 



July 2014 | Volume 5 | Article 371 | 2 



Hu et al 



Ammonia oxidizers and denitrifiers in intertidal wetland 



CO 

II 
C 

d 



o 

"D 
C 

CO 



2! 
a 
ifc 



J2 
O 



"D 
£ 

CO 



0) 
J2 



w o 

.— 

'E 
<u 
■a 



'of 

Z O 
■D 



'of 

Z O 

_i >■ 
■a 



x 



rj> 'o 



O > 



C3> 'o 



Z > 



i . . 

■X O 
If 



-x o 
Z ? 

?! 



'o> o 
II 



X 

Q. 



— IV 

2 I 
in p 









CO 


CO 


CO 


O 


o 


o 


CN 


CN 


CN 




St 




LD 


LD 


LD 


CD 


d 


d 








CO 


CO 


CO 


Csl 


CN 


Csl 














00 




St 


CO 


CO 


CO 


CO 


CO 


Csl 


r^. 






CN 




o 


CO 


CN 


o 


CD 




Csl 


o 






CO 


d 


r-N 




d 




O 


Csi 




CO 


-o 




-H 


+1 


-H 




-H 


-H 


+1 


-H 




-W 


Q 


+1 








o 


CO 




Csl 


o 








CD 




■ 








< ~ > - 


^ 




d 














J 








CO 










CD 




Csl 














•si- 






CD 




CO 




CN 


LD 


CD 


CO 


CO 




Csl 




q 


CN 


CD 




LD 






LD 


CD 




o 


cq 


CO 


CO 


03 


CO 


CO 




csi 


LD 




d 


Csi 




-H 


+1 


-H 


-H 


-H 


-H 


+| 


-H 


-H 


-H 


-H 


+| 


^j- 


CO 


CO 


o 




CD 


LD 


LD 


CD 


CO 


^j- 


CD 


CN 


CO 


p 


p 


*t 










LD 


^; 


CD 




CD 


csj 




CD 


d 




d 




CO 


CO 


CN 


T 


1 


1 












CNJ 


1 


1 


1 




CD 


CD 






CN 












CN 


CO 


CO 




CD 


LD 


CN 


CN 


CD 


LD 


LO 


p 


CD 






csi 




CN 




csi 


Csi 


d 






LO 


CO 


CN 


Csl 


Csl 


Csl 


CN 


Csl 


CN 










CO 


LD 


CD 


*t 


r*» 




LD 


LD 


T— 


CO 


CD 


CO 


p 








CD 


CD 






CD 




LO 












CN 


Csl 


CO 


CO 


o 


CO 


O 


O 


-H 


+1 


-H 


-H 


-H 


-H 


+1 


-H 


-H 


~H 


-H 


+1 


h- 


CD 




Csl 


r-» 


Csl 


CD 




CO 


CO 


Csl 


r^- 


CO 






p 


*t 


LD 


d 


d 


d 


p 


p 












CO 
















<- 






Csl 




<~ 












CN 




O 








LD 




CD 






CD 


o> 






o 


CD 


CN 




CO 


Csl 


o 


LO 


CO 


Csl 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


-fl 




+) 


-H 


-H 


-j-| 


-Lj 






-f| 


-Lj 


-Lj 


























LO 






CD 


CO 


CO 






p 


CD 




CN 


d 


d 


d 


d 


d 


d 


d 


d 




d 






CD 


CD 




o 


^_ 


LD 


CO 


o 


CD 


Csl 


O 


CO 


CO 


Csl 


o 


Csl 


,— 


LD 


o 


,— 


O 


CO 


CO 


CO 


T_ 


d 


r ~ 


d 


o 


d 


d 


o 


d 


d 


d 


d 


-H 


-H 


-H 


-H 


-H 


-H 


-H 


-H 


-H 


~H 


-H 


-H 


CM 


r— 


r-» 


T— 


CD 










r^. 


LO 


t— 




Csl 


LD 




p 


CN 


CN 


CD 


LO 


CD 


p 


CN 


CO 


CN 


CO 


Csj 


CN 


CN 


d 


d 


d 


CN 


Csi 


CO 




CO 




CD 




CD 




CD 




CD 


LO 


ID 


CO 


CO 


CD 






CO 


CO 


Csl 




CN 


LO 


o 




CN 


d 


d 


d 


d 


d 


d 




CN 




CN 


-j-| 


-H 


+) 


-H 


-j-| 


-|-| 


-Lj 


-)-| 




~H 




-Lj 


























p 


p 


LD 




CN 


CD 


CD 


CO 


CD 






LO 


CO 




CO 


CO 




CO 


CO 


CO 


d 


CO 


LO 












*~~ 


















LD 


CO 


*t 


CD 


CD 


CD 


LD 


Csl 


Csl 






p 


p 


LD 


p 






O 




Csl 


p 




O 


o 


O 


O 


O 


O 


o 


o 


o 


o 


o 


o 


o 


-H 


+1 


-H 


-H 


-H 


-H 


+1 


-H 


-H 


-H 


-H 


+1 




CO 


o 


o 


CO 


00 




o 


r^. 


LO 




CO 


LO 


Csl 


o 


CD 


CO 


p 


o 


r-- 




LO 


CD 


p 


CN 


LD 


d 


*t 


r-.' 




CO 


LD 




r>' 


LO 




CO 




*t 


CO 


CD 






CD 


'si" 


CD 




CD 


o 




o 






CO 


o 








O 


O 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


-H 


+1 


-H 


-H 


-H 


-H 


+1 


-H 


-H 


~H 


-H 


+1 


CD 


Csl 


*t 


*t 


CO 


CD 


CD 




Csl 


CD 


LO 


CD 


O 


CN 


CN 


LD 


CO 


O 


CO 




p 


LD 


CD 


00 


CO 






|^ 




CD 


CO 




CO 


CO 






LL 


< 




LL 


< 




LL 


< 


s 


LL 


< 




CO 


CO 




CO 


LO 




CQ 


LO 




CQ 


00 




























CD 


CD 


CD 


a 


C 


EZ 














E 


E 


E 


E 


E 


E 


q3 


CD 


a3 


CD 


CD 


CD 


E 


E 


E 










C 


c 


rir 


rin 


rin 




=3 








=3 












CL 


CO 


CO 


LO 


< 


< 


< 






§ 


LO 


00 


CO 



v- i- -Q 



E 



cz -c -a 



■5 2 



DNA EXTRACTION AND REAL-TIME QUANTITATIVE PCR 

DNA was extracled from 0.25 to 1 g aliquots of each sample in 
triplicate using the UltraClean™ Soil DNA Isolation Kit (Mo-Bio 
Laboratories, USA) as described before (Li et al., 2012). 

Real-time quantitative PCR was performed as described pre- 
viously (Li et al., 2012), with modification. In short, 10 ng of 
DNA extract from each sample was used as a template, and the 
primers and probes are listed in Table 2. The real-time PCR reac- 
tions were carried out by an Mx3000P real-time PCR system 
(Stratagene, USA). The copy number of the 16S rRNA genes of 
the bacteria and Archaea, the ammonia monooxygenase (amoA) 
genes of the AOB and AOA, as well as the nitrous oxide reduc- 
tase ( nosZ) gene of the denitrifiers were determined. To construct 
standard curves, plasmids containing the target gene were quan- 
tified using a NanoDrop 1000 spectrophotometer and were then 
serially diluted in 10-fold steps before qPCR was performed. Each 
reaction was performed in triplicate. The results of each reaction 
are expressed as the copies of each corresponding gene per gram 
dry- weight of soil. 

PCR AMPLIFICATION AND T-RFLP FINGERPRINTING OF THE amoA AND 
nosZ GENES 

The primer set used for amplification of the AOB amoA gene is the 
same as that for qPCR (Table 2). A partial fragment of the AOA 
amoA gene was amplified using the primer set Arch-amoAF and 
Arch-amoAR (Francis et al., 2005), and the nosZ gene was ampli- 
fied using the primer set nosZ2F and nosZ2R (Henry et al., 2006). 
All forward primers were terminally labeled with the reporter dye 
FAM (6-carboxy-fluorescein). The 50u.l reaction mixture con- 
tained 25 ixl of PCR MasterMixture (Takara, Japan), 10 pmol of 
each primer and 20 ng of DNA template. For winter samples, 
no AOA amoA amplicon was obtained. Therefore, a nested PCR 
was employed to obtain sufficient PCR amplicons. The labeled 
amplicons were checked by agarose gel electrophoresis and then 
purified with a Qiagen Gel Extraction Kit (Qiagen, USA). To iden- 
tify the type of enzymes to be used for T-RFLP fingerprinting and 
to determine which enzyme showed the highest resolution and 
least redundancy, an in silico computation was performed before 
digestion. We used Alu\, Hhal, and Taql for the amoA genes of 
AOB, and Afal, Hhal and Msp\ were used for the amoA genes of 
AOA and the nosZ genes of denitrifiers. Each enzyme produced at 
least six T-RFs in all of the samples. In each T-RFLP assay, termi- 
nal fragments (T-RFs) that were not digested with the restriction 
enzymes were considered redundant and were excluded from the 
results matrix. 

Amplicons were digested with restriction endonucleases 
as described in the manufacturer's protocol (Takara, Japan). 
Amplicons of the AOB amoA gene were digested with AZtil and 
Hhal at 37°C for 3h and with Taql at 65°C for 3h, separately. 
Amplicons of the amoA gene of AOA and the nosZ gene of the 
denitrifiers were digested independently with Hhal, Mspl, and 
A/aIat37°Cfor3h. 

After precipitation with ethanol, DNA fragments were size- 
separated using an ABI PRISM 3730 Genetic Analyzer (Applied 
Biosystems). T-RFs from 50 to 600 bp were analyzed with the 
GeneScan software version 3.7 (Applied Biosystems). The rel- 
ative abundance of the individual T-RF was calculated as the 
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Table 2 | The primers, probes and PCR procedures for real-time quantitative PCR used in this study. 



Target gene 



Primer (Probe) Sequence (5 -3 ) 



Thermal profile 3 



References 



Bacteria 16S rRNA 1369F 
1492R 
(1389F) 

Archaea 16S rRNA ARCH1-1369F 
ARCH2-1369F 
1492R 
(1389F) 

AOB amoA amoA-1 F 

amoA-2R 



CGGTGAATACGTTCYCGG 

GGWTACCTTGTTACGACTT 

CTTGTACACACCGCCCGTC 

CGGTGAATACGTCCCTGC 
CGGTGAATATGCCCCTGC 
GGWTACCTTGTTACGACTT 
CTTGTACACACCGCCCGTC 

GGGGTTTCTACTGGTGGT 
CCCCTCKGSAAAGCCTTCTTC 



1 min at 95°C; 40 cycles of 15s at 95°C, Suzuki et al., 2000 
1 min at 56°C 



1 min at 95°C; 40 cycles of 15s at 95°C, 
1 min at 59°C 



1 min at 95°C; 40 cycles of 15 s at 95°C, 
30 s at 57°C, 45 s at 72°C 



Suzuki et al., 2000 



Rotthauwe et al., 1997 



AOA amoA amo-196F GGWGTKCCRGGRACWGCMAC 

amo-277R CRATGAAGTCRTAHGGRTADCC 
(amo-247) CAAACCAWGCWCCYTTKGCDACCCA 



1 min at 95°C; 40 cycles of 15 s at 95°C, 
40 s at 55° C 



Treusch et al., 2005 



nosZ 



nosZ2F 
nosZ2R 



CGCRACGGCAASAAGGTSMSSGT 
CAKRTGCAKSGCRTGGCAGAA 



1 min at 95°C; 6 cycles of 15 s at 95°C, 
30 s at 65°C with a touchdown of 1°C by 
cycle, 30 s at 72°C, 15 s at 80°C; 40 cycles 
of 15s at 95°C, 15s at 60°C, 30s at 72°C 
and 15 s at 80°C 



Henry etal., 2006 



a PCR procedures have been slightly optimized in this study. 



percentage of total peak area in a given T-RFLP profile, and the T- 
RFs that contributed less than 5% were not considered for further 
analysis. 

STATISTICAL ANALYSIS 

All statistical analyses were carried out using the software pack- 
ages vegan and BiodiversityR for the program R (R Development 
Core Team, 2007). Principal component analysis (PCA) was 
employed to explore by T-RFLP fingerprinting the variations in 
the community composition. T-RF profiles of the amoA genes 
of ammonia oxidizers and the nosZ of denitrifiers were ana- 
lyzed independently. To determine the factors that best explained 
the variation in community compositions, redundancy analysis 
(RDA) was performed. The four different sampling seasons and 
three different vegetation types (categorical variables) were com- 
bined to form the variables season and vegetation, respectively. 
The results of all of the tests were examined by the permutation 
test (1000 permutations) at P < 0.05. The dissimilarities among 
all the samples were then calculated using the Bray-Curtis dis- 
tance, and then, non-metric multidimensional scaling (NMDS) 
ordination was used to access the similarities of the microbial 
community composition across different seasons and vegetation 
types. Soil properties and microbial process rates were concate- 
nated to build a secondary matrix. The vector and surface fitting 
of these variables within NMDS ordination were performed, and 
the significance was tested with the permutation test (1000 per- 
mutations). The most significant vectors or surfaces were selected 
to fit against the resulting bi-plots. 

For the purpose of exploring the correlations among com- 
munity size, community composition, environmental variables, 
and microbial process rates, the values for functional gene 
copy numbers, T-RFLP fingerprinting profiles, soil properties 
and nitrification and denitrification rates were transformed into 



dissimilarity matrices using the Bray-Curtis distance measure. 
Correspondences between these matrices were evaluated using 
the Mantel test based on Pearson's product-moment correla- 
tion and 1000 permutations. Before analyses, the qPCR data 
were logio transformed and then divided by each column 
summation to normalize the numerical distribution and pro- 
vide variance homogeneity. Data from summer/autumn and 
from winter/spring were calculated separately during the Mantel 
test. 

RESULTS 

SOIL PROPERTIES AND MICROBIAL PROCESS RATES 

The soil pH values of 12 samples varied between 7.2 and 9.1. 
The salinity changed more dramatically, ranging from 2.54 to 
9.00 mg per g of dry soil (Table 1). In the summer and autumn, 
BF soils had the lowest salinity values compared with the PA and 
SA. In contrast, BF soil salinity was highest in the winter and 
spring. This increase was the result of increased aquatic salinity 
due to the decreased flow rate of the Yangtze River in winter and 
spring. 

The microbial process rates displayed distinct seasonality. The 
net nitrification rates ranged from —17.24 to 121.80 |xg N g _1 
dry soil h~ 1 . All of the samples collected during the summer and 
spring gave negative net nitrification rates (Table 1), which indi- 
cated the denitrifiers had higher activity compared to the nitrifiers 
in those samples. The net nitrification rates in the autumn and 
winter were significantly higher than those in the summer and 
spring (Table 1). In contrast, the potential denitrification rates 
ranged from 0 to 64.09 u,g N g _1 dry soil h _1 . In the summer 
and spring, the values were generally higher than those in other 
samples and showed the reverse trend of the nitrification rates. 
In these two seasons, the potential denitrification rates in BF soil 
were higher than those in the PA and SA soils (Table 1). 
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FUNCTIONAL GENE ABUNDANCE 

The AOB amoA gene abundance ranged from 1.5 x 10 3 to 
3.2 x 10 6 copies g _1 dry soil. Three winter samples displayed the 
lowest values out of the four seasons, which were one to three 
orders of magnitude lower compared with the others (Figure 2). 
However, the level of the total bacterial 16S rRNA gene was sta- 
ble during the study period. In the winter, the copy numbers were 
not at their lowest but were ~ 1 order of magnitude lower than the 
highest level. 

The AOA amoA genes were always higher in copy number than 
the AOB amoA genes in all of the summer and autumn samples. 
Interestingly, in the three winter samples, no AOA amoA gene was 
detected. However, there was no significant decrease in the total 
number of archaeal 16S rRNA gene copies in these samples. To 
confirm this result, three sampling locations were re-sampled in 
January 2010 and the qPCR for the AOA amoA gene was repeated. 
These results confirmed the absence of the AOA amoA in winter. 

In the spring, the AOB amoA abundance was similar to the 
level of the AOA amoA, and the difference was less than one 
order of magnitude. In spring PA, AOB was even higher than the 
AOA. This was the only AOB dominant sample besides the winter 
samples (Figure 2). For both the AOB and AOA amoA gene abun- 
dance, the levels observed in the BF samples were always lowest in 
spring. 

The copy numbers of the nosZ genes ranged from 1.0 x 10 5 
to 9.2 x 10 7 copies g _1 of dry soil. Similar to the AOA and AOB 
amoA, the nosZ genes also showed seasonal specificity; the lowest 
level was in winter. Compared with amoA, the changes in the level 
of the nosZ gene copy number were relatively small. Within the 



same season, the levels of gene copy numbers were almost at the 
same order of magnitude (Figure 2). 

N-CYCLE COMMUNITY COMPOSITION REVEALED BY T-RFLP 
FINGERPRINTING 

For each gene, we retained the result of the enzyme with the max- 
imum T-RFs and the minimum redundancy. These genes were 
Afal for the AOA amoA, Taql for the AOB amoA and Msp\ for 
nosZ. The most prominent T-RFs, which were present in almost 
all the samples, were 303 bp in the AOB amoA, 69 and 172 bp in 
the AOA amoA, and 105 and 126 bp in nosZ (Figure 3). 

The T-RF patterns of these nitrogen-cycling communities 
showed distinct seasonal specificity. For example, the 297-bp T- 
RF of AOB was only present in the summer samples, and the 
186-bp T-RF was only present in the spring samples. Some T- 
RFs were only detectable in samples from the same or similar 
soil type (SA and PA); the 76-bp T-RF of the AOB was present 
in the SA and PA soils but was absent in the BF soils, and the 397- 
bp T-RF of the denitrifiers was present only in the SA soils in all 
four seasons. These seasonal and vegetative specificities were fur- 
ther investigated by PCA and RDA analyses of the T-RFLP data 
(Figure 4). The two PCA axes explained 61.9, 74.2, and 64.7% 
of the community variation for the AOB, AOA, and denitrifiers, 
respectively, and clearly separated the three communities of the 
summer samples from the winter samples. Autumn and spring 
were also separated, but not as clearly. Generally, the autumn sam- 
ples shared more similarity with the summer samples, whereas 
the spring samples were more similar to the winter samples. The 
clearest separation of vegetation type was between BF and the 




BF SA PA BF SA PA 



Soil types 

FIGURE 2 | Average gene copy numbers of 16S rRNA and the functional alterniflora; PA, soil covered by Phragmites australis. Bac 16S, Bacterial 16S 
genes per g of three types of dry soil collected during the four season. rRNA gene; Arch 16S, Archaeal 16S rRNA gene; Bac amoA, amoA gene of 
Error bars represent the SD BF, soil at bare flat; SA, soil covered by Spartina AOB; Arch amoA, amoA gene of AOA; nosZ, nosZ gene of denitrifier. 
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FIGURE 3 | The community composition revealed by terminal 
restriction fragment polymorphism (T-RFLP) analysis of the 
corresponding genes (Bac amoA, amoA gene of AOB; Arch amoA, 
amoA gene of AOA; nosZ, nosZ gene of denitrifier). The colored squares 
represent relative proportions of terminal restriction fragment (T-RF) during 
four seasons in three soil types. BF, soil at bare flat; SA, soil covered by 
Spartina alterniflora; PA, soil covered by Phragmites australis. 



AOB 



two other types of soil (SA and PA). However, there was no obvi- 
ous separation between the SA and PA samples. For the AOB and 
denitrifiers, the main variations could not be associated with dif- 
ferences in the relative abundance of a single or several T-RFs 
(Figure 4). For the AOA, the 123- and 69-bp T-RFs were the 
first two significant factors that determined the main variation of 
the community (Figure 3). The results of the RDA analysis indi- 
cated that all three community compositional variations could 
be strongly explained by the variables of season and soil type 
(Table 3). 

STATISTICAL ANALYSES 

NMDS analysis gave clear separation of the community 
size/composition in the winter/spring and in the sum- 
mer/autumn (Figure 5). Soil sample in the summer and 
autumn shared a greater number of similar environmental 
factors, such as temperature and substrate concentrations, than 
the same samples collected in the winter and spring. To find more 
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FIGURE 4 | Principal component analysis (PCA) of the sampling sites 
based on the composition of the AOB, AOA, and denitrifier 
communities as determined by the T-RFLP analysis. Each vector 
represents an individual T-RF profile. The values in parenthesis indicate the 
percentage of community variation explained by the axes. BF soil at bare 
flat; SA, soil covered by Spartina alterniflora; PA, soil covered by Phragmites 
australis. 



potential relationships between community size/composition 
and environmental factors, we combined two seasons together 
and performed the mantel test separately, and the correlations 
are listed in Tables 4, 5. 
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Table 3 | The conditional effects of environmental variables on the community compositional variation as determined by redundant analysis 
(RDA). 



Environmental 
Variable 




Constrained proportion 3 






F-ratio 




AOB 


AOA 


Denitrifier 


AOB 


AOA 


Denitrifier 


pH 


0.12 


0.06 


0.14 


1.40 (NS) b 


0.64 (NS) 


1.62 (NS) 


Salinity 


0.10 


0.09 


0.08 


1.06 (NS) 


1.04 (NS) 


0.92 (NS) 


NH+ 


0.06 


0.12 


0.03 


0.69 (NS) 


1.33 (NS) 


0.33 (NS) 


NO - 


0 28 


0 21 


0 21 


0.0/ \ I 


? 73 l*\ 


9 79 i*\ 


Total N 


0.10 


0.16 


0.16 


1.06 (NS) 


1.97 (NS) 


1.87 (NS) 


Total C 


0.11 


0.14 


0.11 


1.26 (NS) 


1.67 (NS) 


1.22 (NS) 


C/N ratio 


0.12 


0.07 


0.07 


1.41 (NS) 


0.81 (NS) 


0.81 (NS) 


NiR 


0.13 


0.11 


0.04 


1.55 (NS) 


1.25 (NS) 


0.41 (NS) 


DNiR 


0.19 


0.21 


0.16 


2.30 (*) 


2.68 (*) 


1.93 (NS) 


Season 


0.51 


0.49 


0.44 


2.83 (»»») 


2.54 (*) 


2.10 (*) 


Vegetation type 


0.33 


0.41 


0.46 


2.22 (») 


3.15 (*) 


3.85 (***) 



Abbreviations: NiR, net nitrification rates; DNiR, potential denitrification rates. 

a Constrained proportion represents the proportion of each variable explained by the RDA. 

b NS, not significant at P > 0.05; *0.01 < P< 0.05; ***p< 0.007. 



Our statistical analyses, however, did not find any signifi- 
cant correlation between pH and community size in any season 
(Tables 4, 5). This result was not surprising because the three veg- 
etative types of soils we studied had similar pH values (Table 1). 
Our data showed a significant correlation between the differences 
in the salinity and community size and composition in the sum- 
mer and autumn samples. In the winter and spring samples, only 
the difference in the AOA community composition correlated 
significantly with the differences in the salinity. Other soil prop- 
erties, such as the C/N ratios, significantly correlated with the 
AOB, AOA, and denitrifier community sizes and compositions 
in the summer and autumn. However, similar to salinity, these 
correlations were no longer significant in the winter and spring. 

DISCUSSION 

In many environments, microbial communities display both spa- 
tial and temporal patterns that may be associated with changes 
in the community size and composition. In this study, over four 
seasons, we sampled soil covered by three different types of vege- 
tation (including no vegetative cover) and analyzed which factor 
affects the community most strongly. 

Many previous studies have demonstrated a greater abundance 
of archaeal than bacterial amoA genes in many environments 
(Leininger et al, 2006; He et al, 2007; Lam et al, 2007; Wuchter 
et al., 2007; Nicol et al, 2008; Hong et al, 2013), including estu- 
aries (Beman and Francis, 2006; Abell et al., 2010; Bernhard 
et al., 2010; Li et al., 2012). Our data from the summer and 
autumn samples taken from locations with different vegetative 
cover showed that the AOA were dominant; however, the pattern 
was not the same in the winter and spring samples (Figure 2). 
In the winter, the AOA amoA gene was completely undetected, 
while the AOB amoA gene copy numbers also declined but not as 
dramatically as for the AOA (Figure 2). In March 2007, the AOA 
were also undetectable in the bare flat soil of the Dongtan wet- 
land, the same site used in this study (Li et al., 2012). A slight 
decrease in the AOA abundance in the winter has been previ- 
ously reported in a study of long-term fertilized red soil (He et al., 



2007). The differences in the AOA community size were signifi- 
cantly correlated with differences in the net nitrification rates in 
the summer and autumn (Tables 4, 5). In the spring, although 
the AOA displayed relatively high amoA gene abundance, the net 
nitrification rate still decreased (Table 1 and Figure 2). This result 
suggests that in the winter and spring, the nitrification may be 
principally affected by the AOB instead of the AOA. The seasonal 
variation in the SE USA coastal waters showed that the AOA were 
100 to 1000-fold more abundant during their peak than at other 
times of the year, whereas the abundance of AOB varied < 10- 
fold over the same period (Hollibaugh et al., 2014). Although 
year-round AOA domination in subtropical macrotidal estuary 
sediment in Australia has been recently reported in which the 
AOA community's abundance during the four seasons displayed 
no significant changes (Abell et al., 2010), this may be related 
to the less dramatic seasonal differences in temperature during 
the annual cycle in the Australian system (19-30°C) vs. Shanghai 
(6.0-34. 1°C). A previous study indicated that temperature is one 
of the most important factors controlling the distribution of 
both the AOB and AOA (Sahan and Muyzer, 2008). In combi- 
nation with these previous findings, our results confirmed that 
temperature may be the primary seasonal factor affecting ammo- 
nia oxidizer- community size, and the seasonal impacts on the 
AOA-community size are much greater than the AOB. 

The denitrifier community size, estimated by nosZ gene abun- 
dance, also displayed clear seasonal specificity with the lowest 
level observed in the winter. In the summer and autumn, a signif- 
icant correlation between the potential denitrification rates and 
denitrifier community size was observed. These findings are con- 
sistent with a recent study that showed denitrifier community size 
was related to the corresponding process rate (Patra et al., 2006). 
However, in our study, this correlation was not significant in the 
winter and spring. As was the case with ammonia oxidizers, veg- 
etation had a negligible effect on the denitrifier community size 
within this narrow range of soils. Because the nosZ gene primer 
set could not cover atypical nosZ genes (Orellana et al., 2014) 
and incomplete denitrifiers do not contain the nosZ gene (Henry 
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FIGURE 5 | Non-metric multidimensional scaling (NMDS) ordination of 
sampling sites based on the composition of the AOB, AOA, and 
denitrifier communities as determined by the T-RFLP analysis. The 

seasonal and vegetative patterns of the communities are visualized by 
enclosing the samples from the same season or same soil type in a 
colored ellipse (at 95% confidence). Significant environmental variables 



(P < 0.05) were bi-plotted as vectors. The lengths of the vectors indicate 
the strength of the correlation between the variables and the ordination 
scores. Each vector points to the direction of increase for a given 
variable. Contours represent the net nitrification rate data fit to ordination 
scores. BR soil at bare flat; SA, soil covered by Spartina alterniflora; PA, 
soil covered by Phragmites australis. 



et al., 2006), the denitrifiers detected in this study were only the 
complete denitrifiers that can convert nitrate to nitrogen gas. 

Statistical analysis of the T-RFLP fingerprinting results 
(Figure 5) indicated that seasonal change was the most impor- 
tant indicator for the AOA community composition patterns in 
this wetland ecosystem, and vegetation was more important for 
the AOB community patterns. This finding is in agreement with 
previous studies that suggested vegetation has different effects on 



different microbial communities (Patra et al, 2006; Wang et al., 
2013). The Mantel test results indicate that the functional gene 
abundances significantly correlated with the T-RFLP fingerprint- 
ing results, both in the summer/autumn and in the winter/spring 
(Tables 4, 5), thus implying that there was seasonality in the 
growth or death of the community members. 

Interestingly, in all seasons, the differences in the AOB com- 
munity size and composition were significantly correlated with 
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the differences in the net nitrification rates; and the differences 
in the denitrifier community composition were significantly cor- 
related with the differences in the potential denitrification rates. 
Studies in other soil environments indicated that community 
composition plays a minor role for the process rates when com- 
pared with the role of community size (Hallin et al., 2009). 
However, our results suggest that the community size as well as 
the community composition affects the microbial process rates in 
the environment. 

In this study, distinct changes in the AOB-, AOA-, and 
denitrifier-community sizes and compositions in during all four 
seasons and three vegetative types of soils were found. All of 
the community sizes showed strong seasonal specificities. The 
changes in the community composition showed distinct correla- 
tions based on both the season and vegetation type. Specifically, 
the composition of the AOA community was more significantly 
correlated with seasonal factors and the denitrifier- community 
composition was more significantly correlated with vegetation 
type. A unique absence of AOA in winter was observed, which 
suggests that seasonal factors have the most significant effects on 
this community. Both the community size and community com- 
position contributed to the changes in the process rates, including 
the nitrification and denitrification rates. However, this study 
only included limited sample locations and has been performed 
within a time frame of only 1 year. To understand the issue exten- 
sively, replicate site locations with the analogous vegetative covers 
should be studies over several seasonal cycles. 
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